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INTEGRATED  DETECTOR  TECHNOLOGY  FOR  CORROSION  INSPECTION 


Summary 

This  project  investigated  the  use  of  gas  micro-strip  detectors  as  a  means  of  improving  the 
throughput  of  an  innovative  technique  for  detecting  hidden  corrosion  in  aging  aircraft.  The 
evaluations  included  the  development  of  system,  detector  and  data-readout  concepts  capable  of 
meeting  the  resolution  and  throughput  goals  for  this  instrument,  and  the  modeling  and 
computation  of  x-ray  and  electron  interactions  in  the  detector. 

The  results  of  this  evaluation  indicate  that  these  detectors  can  be  adapted  to  this  application  and 
can  achieve  the  spatial  resolution  and  count  rates  required  to  support  a  system  throughput  up  to 
25  m/hr  for  detection  of  5%  corrosion  in  lapjoints.  This  throughput,  nearly  two  orders-of- 
magnitude  faster  than  the  existing  demonstration  system,  is  sufficient  to  support  the  use  of  this 
instrument  as  a  practical  adjunct  to  existing  technologies.  The  next  step  in  development  would  be 
the  fabrication  of  this  detector  and  the  demonstration  that  the  sensitivity  and  throughput  targets 
can  be  achieved  in  a  demonstration  system. 

Introduction: 

Corrosion  is  the  number  one  maintenance  problem  confronting  the  US  Air  Force.  Corrosion  in 
lap-joints  is  a  problem  of  sufficient  magnitude  to  justify  the  development  of  an  NDE  instrument, 
if  a  reasonably  priced,  high-throughput,  detector  is  available.  This  research  has  been  directed 
toward  developing  and  evaluating  the  technology  required  to  design  and  fabricate  such  a 
detector. 

ARACOR’s  instrumental  concept,  called  COREX,  is  based  on  the  energy-discriminated 
detection  of  x  rays  backscattered  from  the  region  of  suspected  corrosion.  This  research  program 
has  built  upon  a  previous  Air  Force  PRDA  in  which  ARACOR  developed  the  instrumental 
concept  and  upon  an  ongoing  program  through  the  UDRI  to  demonstrate  its  effectiveness  at 
detecting  corrosion  in  aircraft  lap-joints. 

While  the  previous  programs  have  demonstrated  that  the  detective  modality  is  effective  in 
manifesting  low-levels  of  corrosion  products  in  aircraft  lap-joints,  the  demonstration  technology, 
which  uses  discrete  Xenon  proportional  counters  and  NIM  electronics  is  too  slow  and  expensive 
to  be  practical  for  use  as  an  inspection  tool.  In  response,  ARACOR  has  developed  a  concept 
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using  integrated  proportional-counter  array  technologies  developed  for  high-energy  particle 
accelerators. 

Program  Objectives 

Our  goal,  for  this  effort,  has  been  to  evaluate  critical  issues  which  could  lead  to  an  early 
evaluation  of  the  practicability  of  the  integrated-detector  concept  and  which  would  develop  tools 
and  models  to  set  the  strategic  outline  of  the  subsequent  technology  development. 

The  specific  objectives  were: 

1)  Develop  a  model  of  the  x-ray  and  initial  photoeiectron  interactions  in  the  detector: 

For  a  micro-Gap  detector  geometry  used  with  high-pressure  Xenon  gas,  model  absorption 
of  the  incident  x-ray  spectrum,  escape  of  the  Xe-K  fluorescence,  and  primary  photo- 
electron  emission  from  the  substrate.  Determine  the  pulse-height  spectrum  and  the 
measurement  throughput. 

2)  Develop  a  model  of  the  electron  avalanche  and  ion  transport  in  the  detector: 

For  a  given  geometry  {i.e.  size  and  spacing  of  the  electrode  structures)  and  gas  pressure, 
calculate  the  gain  and  temporal  response  of  a  unit  cell  of  the  detector. 

Technical  Results 

Since  the  original  proposal  for  this  research  was  written,  the  technological  context  of  this  work 
has  continued  to  evolve.  In  particular,  one  technology,  eddy  current,  has  emerged  as  the  Air 
Force’s  primary  instrument  for  detection  of  lap-joint  corrosion  in  aging  aircraft.  This  does  not 
obviate  the  need  for  the  COREX  technology;  the  “false-positive”  indication  rate  of  the  eddy- 
current  instruments  is  high  enough  that  there  is  a  pressing  need  to  verify  the  presence  of 
corrosion  before  expensive  repairs  are  undertaken.  However,  as  a  secondary  instrument,  the 
emphasis  for  the  COREX  instrument  will  be  more  on  keeping  production  costs  low  than  on 
achieving  high  inspection  throughput.  While  this  has  colored  our  thinking  about  what  this 
instrument  will  look  like,  the  major  cost/throughput  tradeoff  decisions  will  take  place  in  the 
development  of  the  prototype;  the  proposed  use  of  an  integrated  detector  technology  still  appears 
to  be  the  best  way  to  minimize  fabrication  costs. 

The  detector  technologies  discussed  in  the  proposal,  the  micro-strip  gas  chamber  (MSGC),  the 
micro-gap  chamber  (MGC),  and  the  various  mesh-based  techniques  have  continued  their 
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development  and  evaluation. Of  these,  the  MGC  appears  to  be  best  for  our  application  due  to 
its  greater  fabrication  simplicity,  freedom  from  charging  effects  and  physical  robustness.  We 
have  used  this  technology  as  our  baseline,  however,  we  are  spectators  in  the  rapid  development 
of  these  technologies.  Any  of  these  could  most  likely  be  adapted  for  use  with  our  concept  with 
some  modification  to  our  design  and  some  customization  of  the  technology  “as  developed.” 


Drift  Cathode  Pla 


Anode  Cathode 


Microstrip  Gas  Chamber 


MicroGap  Gas  Chamber 


Figured.  The  elements  of  the  MSGC  and  MGC  are  shown  here  {not  to  scale).  The  space 
between  the  electrode  structure  and  the  Drift  Cathode  Plane  is  filled  with  an 
absorbing  counter  gas  (typically  Ne,  Ar,  or  Xe  mixed  with  a  quench  gas).  Typical 
spacing  (Anode  to  Anode)  is  200  pm,  while  the  drift  gap  is  typically  3-10  mm. 


The  electrode  structures  of  the  MSGC  and  MGC  are  illustrated  in  Figure  1.  In  both  detectors,  the 
proportional-counter  avalanche  occurs  in  the  small  (A-K)  gap  between  the  anode  and  cathode. 
Because  minimal  ion  current  flows  back  to  the  drift  cathode,  most  of  the  current  is  captured  very 
quickly.  One  problem  with  the  MSGC  has  been  the  charging  of  the  substrate  in  the  A-K  gap.  The 
design  of  the  MGC  largely  obviates  this  problem. 

The  strategic  approach  for  this  program  began  by  laying  out  an  optimal  system  geometry  and  a 
definition  of  the  x-ray  and  detector  optics.  We  then  built  an  x-ray  model  to  evaluate  the 
performance.  The  desired  performance  characteristics,  and  the  detector  geometry  raised  a 
number  of  feasibility  questions,  particularly  the  feasibility  of  simultaneously  meeting  the  energy 
and  spatial  resolution  without  giving  up  the  open  geometry  and  count  rate.  In  addressing  the 
feasibility  questions,  alternatives  and  modifications  were  introduced  until  we  achieved  an 
optimized  solution. 
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System  and  Detector  Configurations 


The  final  configuration  of  the  system  is  shown  to  scale  in  Figure  2.  A  number  of  issues  for  the 
detector  design — collimation  angles,  drift  distances,  spatial  resolution,  gas-pressure,  x-ray 
capture  and  escape — needed  to  be  evaluated  in  a  specific  geometrical  context  to  be  meaningful. 
The  x-ray  source  is  based  on  that  which  has  been  used  in  the  previous  research.  To  maximize  the 
throughput,  a  compact  geometry  was  used  and  the  acceptance  angles  of  the  detector  were  opened 
up  a  much  as  possible  within  constraints  posed  by  the  requirements  for  spatial  and  angular 
resolution.  It  is  expected  that  a  second  detector  package  will  be  located  to  the  left  side  of  the 
source.  It  was  left  out  of  the  figure  for  graphical  clarity. 


Figure  2.  This  cross-section  of  the  system  shows  the  spatial  relationships  between  the 
components  and  the  optical  constraints  of  the  collimators.  The  system  scans  a  swath 
along  the  lap-joint  parallel  to  the  plane  of  this  view.  The  x  rays  are  collimated  in  a 
fan  which  spreads  across  the  lap  Joint;  the  illuminated  surface  is  a  4-mm  wide  strip 
which  runs  normal  to  the  page.  X  rays  scattered  from  the  lap-joint  volume  are 
collimated  into  planes  parallel  to  the  page.  The  dashed  lines  show  the  umbral  and 
penumbral  extents  of  the  collimated  beams. 

The  performance  goals  for  the  system,  and,  by  derivation,  the  detector,  are  based  on  experience 
with  a  two-channel  breadboard  system  which  has  been  fabricated  and  tested  at  ARACOR.  The 
performance  objectives  for  the  new  system  are: 

Detection  Standard:  The  system  is  evaluated  in  terms  of  its  ability  to  detect  the  oxidation  of  5% 
of  the  aluminum  in  a  lap-joint  made  from  two  1-mm  thick  sheets. 
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Spatial  Resolution:  The  existing  breadboard  system  views  a  volume  element  approximately  4- 
mm  in  diameter.  The  images  sample  the  inspection  region  on  a  2-mm  grid;  this  appears  to 
be  about  optimal.  The  goal  was  for  the  new  system  to  match  these  values. 

Throughput:  Using  two  detectors,  the  breadboard  system  takes  approximately  1/2  second  to 
obtain  a  single  data  point  with  a  signal-to-noise  ratio  (SNR)  of  1:1  for  the  detection 
standard.  During  the  evaluations  referred  to  above,  we  found  that,  at  an  SNR  of  1:1, 
patterns  of  corrosion  could  be  confidently  identified.  The  goal  was  to  increase  the 
throughput  by  1-2  orders  of  magnitude. 


The  detector  concept  is  illustrated  in  Figure  3.  The  parameters  in  this  concept,  (materials, 
potentials,  and  dimensions)  have  been  chosen  to  meet  the  performance  goals  of  the  system. 


Active  Detector  Volume  Drift  Cathode  (-2500  V) 


Dead  Volume 
Resistive  Field-Plate  (Si) 


Figure  3.  The  detector  concept  is  shown  from  the  same  direction  used  in  Figure  2.  In  this  design, 
no  radiation  impinges  on  the  cathode  or  MGC  surfaces — greatly  reducing 
photoemission  problems.  The  resistive  field  plates  keep  the  electric  potential  surfaces 
in  a  plane-parallel  geometry.  The  active  volume  is  similarly  constrained  by  resistive 
field  plates  along  the  surfaces  parallel  to  the  page. 


The  X-ray  Performance  Model 


The  signal  received  in  the  detector,  for  x  rays  of  coherently-scattered  from  aluminum  in  the 
volume  element  (voxel)  being  inspected,  can  be  written. 


^Al,coh 


=  Le 


^vox  Al.coh^^scut^ 

D]  do. 


(1) 
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where  explicit  energy  dependence  has  not  been  shown.  The  terms  in  this  equation  are: 


K 

keV 

The  energy  of  the  x  ray  as  it  leaves  the  source  and  is  incident 
on  the  sample. 

^out 

keV 

The  energy  of  the  x  ray  as  it  leaves  the  sample  and  is  incident 
at  the  detector. 

^Al,coh(^in) 

X  ray/keV-sr-s 

Detected  signal. 

lo 

X  ray/keV-sr-s 

X-ray  source  spectrum. 

Min(Ein) 

cm^/g 

Absorption  coefficient  for  x  rays  entering  the  volume  element 
and  being  attenuated  by  the  overlying  sheet  of  aluminum. 

Mout(E„„,) 

cm^/g 

Absorption  coefficient  for  x  rays  leaving  the  volume  element 
and  being  attenuated  by  the  overlying  sheet  of  aluminum. 

Pt 

g/cm^ 

Product  of  density  and  thickness  for  the  overlaying  layer. 

rad 

Angle  of  entrance  to  the  voxel  (0  for  this  model). 

9out 

rad 

Angle  of  exit  from  the  voxel  (37t/4  for  this  model). 

rad 

Scattering  angle  (37i/4  for  this  model). 

nivox 

g 

Mass  of  aluminum  in  the  voxel. 

cm 

Distance  from  the  source  spot  to  the  voxel, 

^(^Ai.coh(Em,6 

i  ) 

Partial  differential  scattering  cross-section  for  coherent  scatter 

from  aluminum. 

do. 

-  cm  /g  oi 

sr 

Solid  angle  of  the  detector  as  seen  from  the  voxel,  and 

^det(^out) 

Quantum  efficiency  of  the  detector. 

Angular  variations  in  the  scatter  of  the  x  rays  must  be  tracked  accurately;  the  scattering 
distributions  in  this  energy  range  show  very  strong  angular  dependence,  with  the  coherent  (also 
called  Rayleigh)  scatter  being  strongly  peaked  in  the  forward  direction.  The  signal  detected  for 
incoherent  (also  called  Compton)  scatter  and  for  scatter  from  oxygen  can  be  similarly  written 
with  a  change  in  the  choice  in  partial  differential  cross  section.  For  coherent  scatter,  the  x-ray 
energy  is  unchanged  by  the  scattering  process.  However,  there  is  an  energy  loss  in  incoherent 
scattering. 


6 


(2) 


£  _ Zh _ 

where  is  the  rest  mass  of  the  electron  in  keV. 

The  X  rays  interact  in  the  detector  either  through  coherent  scatter,  Compton  scatter  or 
photoelectric  absorption.  The  coherent  scatter  is  elastic;  no  energy  is  deposited.  Compton  scatter 
deposits  a  small,  angularly  dependent,  energy  in  the  detector.  The  energy  deposited  is  below  the 
discrimination  threshold  and  is  not  tracked  by  the  program.  It  will,  however,  contribute  a  small 
low-energy  tail  to  the  spectrum.  The  main  interaction  in  this  energy  region  is  through  the 
photoelectric  effect  which  ejects  an  electron  from  an  atom  of  the  gas.  The  energy  of  the 
interacting  x  ray  is  split  between  the  binding  energy  of  the  electron  and  the  kinetic  energy  of  the 
ejected  electron;  the  electron  deposits  its  kinetic  energy  in  the  gas.  The  hole  left  in  the  inner  shell 
(mostly  either  K  or  L)  of  the  gas  atom  will  decay,  mostly  through  K-  or  L-shell  fluorescence  or 
by  emitting  an  Auger  electron.  Inter-shell  energy-transfers  play  a  minor  role  in  the  detector 
response.  The  energy  of  an  Auger  electron  is  quickly  absorbed  by  the  counter  gas.  Fluorescent 
X  rays  may  be  captured  or  may  escape  from  the  counter.  Only  single  x-ray  interactions  (and  the 
subsequent  decays)  are  tracked  by  the  model. 

The  flow-logic  of  the  x-ray  performance  model  is  illustrated  in  Figure  4.  The  software  is 
implemented  in  Mathematica®  3.0  with  the  front  end  running  on  a  Macintosh  G3  PowerPC  and 
the  kernel  running  on  an  SGI  Indigo2  workstation.  In  operation,  the  program  is  set  up  with  an 
appropriate  set  of  parameters  and  run.  A  number  of  intermediate  calculations  and  data  plots  are 
output  during  the  run,  culminating  in  a  calculation  of  the  time  required  for  a  single  detector 
channel  to  obtain  an  SNR  of  1:1  for  the  detection  standard.  Since  the  total  running  time  for  the 
program  is  about  80  s,  it  is  not  difficult  to  evaluate  the  effects  of  varying  any  particular 
parameter. 

The  x-ray  libraries  used  contain  a  number  of  x-ray  related  data  bases  and  functions  which  have 
been  accumulated  by  the  Principal  Investigator  over  a  number  of  years  and  implemented  as 
Mathematica-callable  functions.  The  various  x-ray  cross-sections  are  derived  from  the  Lawrence 
Livermore  National  Laboratories  EPIC  x-ray  data  tables. 

To  keep  the  evaluation  time  down,  the  spectrum  is  broken  into  discrete  (0.5  keV)  bins.  All 
interaction  coefficients  and  functions  are  then  defined  at  the  bin  energies;  this  avoids  repetitive 
interpolations  from  the  data  base.  This  also  presents  a  complication  for  the  calculation  logistics. 
In  general,  x  rays  which  scatter  incoherently  do  not  lose  energy  in  multiples  of  0.5  keV.  In 
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addition,  the  width  of  the  energy  bins  changes  after  the  scatter.  The  step  in  Figure  4  called 
Generate  “Pre-scatter"  Spectrum  divides  the  source  into  unequal,  and  unequally-spaced,  bins  so 
that  the  x-ray  energies  will  exactly  fit  into  the  final  bins  after  the  Compton  scatter. 


Input  X-ray  Libraries 
Define  Scattering  Functions 

i 

Define  Detector  Model 

i 

Input  System  Parameters 
Define  Corrosion  Model 

y 

Generate  Source  Spectrum  Generate  "Pre-scatter"  Spectrum 

I  i 

Apply  Coherent  Scatter  Apply  Incoherent  Scatter 

\  y 

Merge  Scattered  Spectra 

i 

Apply  Detector  Model 

i 

Evaluate  Signature 

Figure  4.  The  calculational  flow  of  the  detector-physics  model 


The  four  scattering  components  (coherent  and  incoherent  from  aluminum  and  oxygen)  are 
computed  separately  and  then  merged  into  a  final  “scattered”  spectrum.  Its  interaction  with  the 
detector  physics  is  then  evaluated.  At  each  point  in  the  sample  a  “coherent”  signal  S,.  and  an 
“incoherent”  signal  Sj  are  measured.  If  we  define  these  signals  (counts-per-second)  for 
uncorroded  aluminum  (S^.,,  S^i.,)  and  for  the  oxidized  aluminum  (S^q.,,  S^q.,)  then  we  can 
measure  the  so-called  Rayleigh/Compton  (or  RC)  ratio  for  the  two  samples. 

=  (3) 

RC  =  ^Aio-c/  r4\ 

A/o-,' 


The  instrument  operates  by  scanning  a  region  of  a  lap-joint;  the  data  are  presented  as  an  image  of 
the  RC  ratio.  The  presence  of  corrosion  is  manifested  by  a  pattern  of  dark  splotches  on  the  image 


8 


of  the  lap-joint  region.  The  signature  of  corrosion  can  thus  be  expressed  as  a  ratio  of  the  signals 
in  the  corroded  and  uncorroded  regions. 


Sig  = 


(5) 


Because  the  contrast  of  an  inspection  image  on  a  monitor  can  be  arbitrarily  enhanced,  the 
limiting  factor  in  the  sensitivity  is  set  by  the  ratio  between  the  signature  and  the  image  noise 
which  is  determined  by  the  Poisson  statistics  of  the  signals.  For  a  sample  time  t^.,  the  standard 
deviation  of  the  RC  ratio  (for  a  single  detector)  can  be  found  from: 


G 


RC  ~ 


1 


(6) 


As  discussed  above,  when  the  SNR  becomes  1:1  (Sig=aRc).  the  presence  of  corrosion  becomes 
evident  in  the  image.  For  two  detectors,  and  for  the  final  system  configuration,  the  sample  time 
to  detect  5%  corrosion  is  found  to  be: 


2g 


■  +  ■ 


=  0.28j 


RC 


(7) 


This  suggests  that  a  system  designed  according  to  this  model  could  scan  at  25  m/hr  (sampling  on 
a  2-mm  grid)  and  detect  5%  corrosion.  In  comparison  with  the  existing  2-detector  breadboard, 
the  sample  time  is  only  about  20%  shorter.  However,  because  of  its  integrated  detector  design, 
the  detector  simultaneously  measures  40  points  across  an  80-mm  lap-joint,  effecting  a  speedup 
factor  of  50!  This  completely  meets  our  design  goals. 

The  performance  goals  for  the  system  have  been  met  by  significantly  opening  up  the  collimation 
on  the  system.  In  particular,  the  distance  between  the  drift  cathode  and  the  MGC  electrode 
structure  is  3-4  times  larger  than  existing  detector  designs  (see  Figure  3).  In  our  original 
proposal  we  were  concerned  with  the  effects  of  x  rays  directly  impinging  on  the  drift  cathode 
and  MGC  electrodes.  We  largely  eliminated  this  problem  by  bringing  in  the  x  rays  parallel  to  the 
electrodes,  so  that  the  direct  irradiation  was  avoided.  In  the  next  sections,  we  will  look  at  the 
implications  of  these  decisions. 


The  electron  Cloud 


As  the  X  rays  are  absorbed  in  the  counter  gas  they  generate  fast  electrons  which  ionize  more  gas 
molecules  as  they  scatter  and  slow  down.  To  measure  the  energy  of  the  interacting  x  ray  it  is 
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necessary  to  collect  the  complete  cloud  of  electrons,  thus  we  need  to  know  how  big  it  is.  The  RC 
signature  is  determined  largely  by  the  incoherent  scatter  peak  near  30  keV;  this  generates 
photoelectrons  in  the  gas  principally  in  the  15-25  keV  range.  To  determine  the  size  of  the 
generated  cloud,  a  monte-carlo  simulation  of  the  electron  trajectories  was  used  to  define  the 
initial  size  of  the  charge  cloud  created  by  the  photoelectrons.  This  program,  NIST  Micro  MC, 
was  developed  as  a  software  tool  for  the  calculation  of  the  characteristic  x-ray  emission  under 
electron  bombardment  and  is  used  to  model  problems  in  electron  probe  x-ray  microanalysis  and 
analytical  electron  microscopy. An  image  of  a  simulation  of  500  electrons  interacting  with  the 
gas  is  shown  in  Figure  5.  In  this  figure,  the  electrons  enter  at  the  red  mark  (top  center).  We  can 
see  that  few  electrons,  and  very  little  energy,  is  deposited  more  than  about  0.75  mm  from  the 
point  of  generation  of  the  photoelectrons. 


Eo(kV)=25;  Specimen=Xe;  TlltIdeg]=0:  Traj.No.=500 


1 000  microns  BS  coefficient  =  AB.AWi 


Figure  5.  A  simulation  of  25  keV  electrons  entering  the  counter  gas  (2.5  bar)  and  scattering.  As 
the  electron  loses  energy  the  color  shifts  from  red  to  black. 

As  each  electron  track  is  accelerated  to  the  MGC  electrode  structure  it  will  continue  to  diffuse. 
To  evaluate  the  magnitude  of  the  diffusion,  we  used  results  from  the  Magboltz  program  which  is 
available  from  the  CERN  web-site  at  http://consult.cem.ch/writeup/garfield.  This  program  does 
Monte  Carlo  simulations  of  charge  drift  and  diffusion  for  complex  mixtures  of  counter  gases  in 
support  of  gas-counter  design  for  the  high-energy  physics  community.  From  this  software  we 
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were  able  to  obtain  diffusion  coefficients  for  drift  parallel  to  [Dp=  0.52  (cm/Bar)*^^]  and 
transverse  to  [Dt=  0.42  (cm/Bar)'^]  the  electric  field  for  a  mixture  of  90%  Xe/10%  CO2. 
Broadening  of  the  distribution  can  be  found  from: 

FWHM^  =  135a  ^  =  2350^“"^/  =  0.53mm  (8) 

FWHM,  =  235a,  =  2.35D,  =  0.43mm  (9) 

To  estimate  the  effect  of  diffusion  on  the  breadth  of  the  charge  track,  these  diffusion  widths  can 
be  convolved  with  the  widths  of  the  charge  tracks;  this  will  lead  to  only  a  slight  broadening  of 
the  distribution.  In  general,  we  find  that  the  charge  tracks  will  be  up  to  about  0.9-nim  in  both 
directions.  This  has  implications  for  the  temporal  response  of  the  system  and  for  the  charge- 
collection  geometry.  The  temporal  response  is  limited  by  the  time  it  takes  to  collect  the  full 
charge  cloud.  From  the  Magboltz  program  were  also  able  to  obtain  the  drift  velocity  (V(yjj=0.45 
cm/ps)  for  electrons  in  the  0.30  V/cm-torr  electric  field.  The  collection  time  t^„„  is  then: 

=0.21X5  (10) 

The  peaking  time  of  the  pulse  amplifier  should  be  no  shorter  than  the  collection  time. 

In  order  to  correctly  measure  the  energy  of  the  x  ray,  it  is  necessary  to  collect  all  of  the  electrons 
resulting  from  the  x  ray’s  electron  cloud.  Our  strategy  for  making  this  measurement  is  illustrated 
in  Figure  6.  Groups  of  MGC  anodes  are  connected  together  so  that  they  collect  the  proportional- 
counter- amplified  signals  over  a  1-mm  wide  strip.  After  these  signals  are  amplified  and  shaped 
by  a  fast  charge  preamplifier/pulse-amplifier  combination,  they  are  summed  with  each  adjacent 
strip.  If  a  given  adjacent-strip  sum  is  greater  than  its  neighbors  and  fits  within  one  of  the  two 
discriminator  bands  used  to  select  either  the  coherent  or  incoherent  signal,  a  count  is  registered  in 
the  appropriate  register. 

The  problem  addressed  by  the  data-readout  strategy  is  that  of  partially-collected  charge  clouds. 
X  rays  in  the  critical  25-30  keV  range  will,  in  general,  deposit  their  signal  in  two  adjacent  strips. 
We  know  from  the  charge-diffusion  studies  that  it  will  be  very  rare  for  these  charges  to  be 
distributed  into  more  than  two  strips.  Thus,  the  comparison  of  adjacent  strip  sums  will  select  the 
pair  which  captured  the  full  charge.  It  will  then  be  classified  appropriately. 
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Diffused  Drift  Pattern  from 


MGC  Anodes  (multiplexed)  Pulse  Processing  Decode  Logic  System  Interface 


Figure  6.  The  data  readout  logic  is  designed  to  collect  the  signal  from  a  diffused  charge  cloud 
and  determine  whether  it  fits  within  one  of  the  defined  scatter  categories. 

The  charge-collection  strategy  will  impact  the  energy  resolution  of  the  instrument.  Based  on  the 
Magboltz  results,  we  can  expect  the  performance  of  this  component  of  the  energy-resolution  to 
be  very  good  below  about  30  keV.  At  higher  energies,  the  charge  clouds  become  larger  and,  for 
some  interactions,  a  small  part  of  the  signal  will  not  be  collected.  This  loss  of  resolution  at  higher 
energies  will  not  significantly  affect  the  instrument  performance  as  the  higher  energy  x  rays  are 
not  used  in  the  imaging. 

The  energy  resolution  is  also  affected  by  the  properties  of  the  counter  gas.  It  has  been  found  that 
the  MGC  structures  have  about  20%  better  energy  resolution  than  standard  proportional  counters 
and  multi -wire  proportional  counters  (MWPCs).’  The  small  increase  in  counter-gas  pressure 
from  2  to  2.5  Bar  (that  is  in  comparison  with  the  breadboard  instrument)  will  tend  to  offset  this 
improvement  as  shown  in  pressure  vs.  resolution  evaluations  done  for  hole-structures.  For  our 
application,  good  energy  resolution  is  useful  in  maintaining  image  contrast,  but  not  critical. 

Electronics  Considerations 

The  grouping  of  the  anodes  increases  the  input  capacitance  of  the  detector  to  about  100  pF'  and 
presents  a  potential  issue  for  the  ability  of  the  electronics  to  read  out  the  signal  with  an 
acceptable  noise  level.  To  verify  that  our  readout-noise  goals  were  reasonable  we  used  a 
published  model  for  the  noise-performance  of  particle-detector  electronics.  As  inputs  to  the 
model,  we  used  the  following  parameters: 

Peaking  time  (Semi-Gaussian  shaping);  0.2  ps. 

Total  input  capacitance;  150  pF, 
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FET  transconductance:  17  mS, 

Total  leakage  current:  10  pA, 

Total  parallel  resistance:  100  MQ, 

Temperature:  300  K. 

For  this  type  of  detector,  and  at  a  relatively  fast  peaking  time,  the  performance  is  essentially 
driven  by  the  input  capacitance  and  the  transconductance  of  the  FET.  For  this  model,  we 
estimated  a  readout  noise  of  1665  electrons.  Assuming  a  targeted  avalanche  gain  of  1000,  this 
noise  level  is  equivalent  to  approximately  36  eV  referred  to  the  x-ray  energy.  This  is  negligible 
in  comparison  to  the  -1000  eV  resolution  required. 

Photoemission  Effects 

The  original  concept  for  the  detector  had  the  x  rays  entering  through  the  drift  cathode  and  exiting 
through  the  MGC  electrode  structure.  The  consequences  of  photoemission  from  either  electrode 
surface  were  of  some  concern.  In  the  revised  concept,  the  x  rays  enter  through  the  silicon 
resistive  field  plates.  We  modeled  the  photoemission  from  the  entrance  surface  to  evaluate  the 
contributions  to  the  detector  spectrum.  The  models  extrapolated  primary  (>  50  eV)  and 
secondary  (<  50  eV)  photoemission  data  which  had  been  taken  for  x-ray  energies  from  2.6  to  5.4 
keV.*’  Lower-energy  data  were  not  used  because  their  primary  interaction  was  with  Si-L  shell 
which  would  not  support  extrapolation  to  higher  energies. 

Extrapolation  of  the  primary  photoemission  followed  the  form,'® 

fp<<Ex)\x^^{Ex)  (11) 

where  the  function  fp(Ex)  is  a  slowly  varying  function  of  the  x-ray  energy  which  approximately 
follows  the  energy-dependent  range  of  the  most  energetic  photoelectrons  produced  by  the  x 
rays.  The  function  was  empirically  fit  to  the  form 

fp{Ex)  =  A  +  B*Ex-¥C*Ex^  (12) 

The  bulk  of  the  primary  photoemission  is  emitted  after  scattering  several  times  within  the  solid, 
leading  to  a  fairly  uniform  distribution  of  energies  below  the  “no-loss”  escape  peak. 

Extrapolation  of  the  secondary  photoemission  followed  the  form, '® 

(13) 
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where  the  function  fj(Ex)  is  a  very  slowly  varying  function  of  the  x-ray  energy  which  is  more 
strongly  weighted  toward  the  lower-energy  photoelectrons  produced  by  the  x  rays.  The  function 
was  empirically  fit  to  the  form: 

f{Ex)^D*Ex^*  (14) 

The  secondary  photoemission  is  emitted  as  a  low-energy  peak  (~l-2  eV).  The  results  of  the 
photoemission  model  are  shown  in  Figure  7.  While  the  model  is  probably  only  accurate  within  a 
factor  of  a  few,  the  contribution  (145  electrons/s)  is  a  factor  of  1000  lower  than  the  signal  due  to 
x-ray  interactions  in  the  detector.  The  photoemission  calculation  is  included  as  Appendix  1. 


Figure  7.  The  modeled  photoemission  from  the  silicon  field-plate  at  the  entrance  to  the  detector 
is  shown.  The  vertical  axis  is  electrons/bin,  the  horizontal  axis  is  the  electron  energy 
in  keV.  Contributions  to  the  overall  detector  spectrum  are  completely  negligible. 

Contributors 

The  majority  of  the  technical  work  and  models  were  done  by  the  Principal  Investigator,  Dr.  Jerel 
Smith.  The  monte-carlo  calculations  of  the  electron  interaction  in  the  gas  were  provided  by  Mr. 
Edward  Franco.  The  evaluations  of  the  electronic  noise  were  contributed  by  Dr.  June  Wang. 

Conclusions 

This  project  investigated  the  use  of  integrated  proportional-counter  array  technologies  to  achieve 
rapid  detection  of  corrosion  in  aging  aircraft  lap-joints  using  ARACOR’s  COREX  inspection 
technology.  The  detector  technology,  developed  for  use  in  high-energy  physics  tracking 
detectors,  promised  good  spatial  resolution  and  high-speed  energy  discrimination  based  on  the 
state  of  the  art.  Extending  this  technique  to  the  detection  of  corrosion  requires  that  the  detectors 
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operate  at  higher  x-ray  energies,  with  larger  inter-electrode  spacings,  and  at  reduced  spatial 
resolution.  In  this  project  we  developed  models  for  the  x-ray  interactions  in  the  detector  and  the 
subsequent  electron  transport  that  govern  the  spatial  resolution  and  temporal  response  of  the 
detector.  This  program  has  been  extremely  successful,  not  only  in  accomplishing  the  formal 
goals  of  developing  physics  models  and  evaluating  the  potential  for  this  technology,  but,  more 
importantly,  in  developing  a  detector  concept  that  can  support  nearly  two  orders-of-magnitude 
increase  in  the  inspection  speed. 


This  project  developed  an  initial  design  for  a  gas  micro-strip  detector  that  achieves  the 
performance  objectives  of  the  COREX.  We  demonstrated  analytically  that  this  detector, 
operating  at  an  inter-detector  spacing  of  2  mm,  can  achieve  count  rates  above  100,000  Hz. 
Further,  we  estimate  that  the  energy  resolution  of  1  keV  is  more  than  adequate  for  resolving  the 
characteristic  signals  from  coherent  and  incoherent  scatter.  The  basic  conclusion  of  this  project  is 
that  this  detector  technology  can  provide  the  basis  for  a  corrosion  detection  system  with  a 
throughput  up  to  25  m/hr.  for  the  detection  of  5%  corrosion  in  lap-joints.  Detection  of  corrosion 
down  to  1%  can  be  accomplished  at  lower  scan  speeds. 

Following  this  successful  demonstration  of  feasibility,  we  recommend  that  the  micro-strip 
detector  be  fabricated  and  experimentally  evaluated  in  a  breadboard  corrosion-monitor 
instrument.  This  will  provide  direct  evidence  of  the  projected  sensitivity  and  system  throughput. 
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■  Appendix  1 :  AFOSR  Detector:  Photoemission  Evaluation 

The  objective  is  to  estimate  the  effect  of  photoemission  from  internal  detector  surfaces  on  the  low-energy  end  of  the  spectral 
response.  In  this  design,  the  x  rays  enter  from  a  Be  window  followed  by  a  0.1  mm  Si  sheet,  and  end  in  a  0.1  mm  Si  sheet  after 
passing  through  5  cm  of  2-atm  Xe  gas. 

<<  "xSec.m" 

xSec.m  Version  1.1.3,  ®  1998  Jerel  A.  Smith 

■  Define,  Eout,  the  spectrum  into  the  detector  ( Specin),  and  the  detected  spectrum  (SpecDet) 

(taken  from  ’^Detector  60  kV,  2.5  bar,  2-mm") 

Ll  =  Li  stPlot  [Transpose  [  {Eout,  Specin}],  PlotRwge ->  {{0,  MaxEi),  (0,  [  Specin]  }}  , 

DisplayFunction  ->  Identity]  ; 

L2  s  ListPlot  [Transpose  [  {Eout,  SpecDet}],  PlotR2mge ->  {{0,  HaxEi},  {0,  Meoc  [Specin]  }}  , 

Plot  Joined  ->  True,  DisplayFunction  ->  Identity]  ; 

Show[{Ll,  L2},  DisplayFunction ->  $DisplayFunction]  ; 


The  primary  photoemission  from  the  silicon  field  plate  can  be  estimated  from  the  smith/bemstein  measurements.  Assume  that 
the  primary  emission  is  smeared  uniformly  from  the  exciting  photon  energy  to  the  lowest  energy  bin.  To  extrapolate  the 
primary  (Yp)  and  total  (Yt)  yields ,  use  an  empirical  power  law  for  the  electron  ranges: 

ExTest =  {2.64,  4.12,  4.55,  4.95,  5.41}  ;  (*  start  edDove  the  k-edge*) 

YpData=  {5.8,  3.5,  3.15,  2.9,  2.63}  ★  10’^  ; 

YxEMuData=  {2.1,  1.85,  1.78,  1.75,  1.82}  *10'®  ;  (*  divided  by  Ex  Mu  (Ex)  *) 

■  Scale  the  Yp  data  with  the  energy-absorption  and  a  power  law  in  Ex 


YpMu  =  YpData  /  xmiAbs  [Si,  ExTest] 

{4.09533  x  10'*,  8.35144  x  10'*,  9.89986x  10'*,  0.0000115531,  0.0000134353} 
Lpd  =  ListPlot  [Transpose  [  {ExTest,  YpMu}],  DisplayFunction ->  Identity]  ; 
pRange[x_]  =  Fit  [Treuispose  [  {ExTest ,  Y^Mu}],  {1,  x,  x^},  x] 

6.47859  X  10'''  +  2.90772  x  lO'"^  x  +  3 . 84026  x  lO*''  x^ 
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LpFit  =  ListPlot [Transpose [{ExTest,  pRango [ExTest] }] /  Plot Joined  ->  True, 
DisplayFunction  ->  Identity]  ; 

Show[{Lpd,  LpFit},  DisplayFunction ->  $DisplayFunct ion]  ; 


Yp[Ex_]  :=  pRauige  [Ex]  iituAbs[Si,  Ex]  ; 

Yp [ExTest]  *  10^  (★  check  the  fit  to  the  original  data  *) 

{5.79529,  3.50546,  3.15678,  2.88586,  2.63496) 

■  This  is  the  **unsmeared”  primary  photoemission 

Sp  =  Table  [Spec  In  I  ij  Yp  [Eout  |[i]]  ]  ,  {i,  nPts}]  ; 

ListPlot [ 

Transpose  [  {Eout,  Sp}],  PlotRange ->  {{0,  MaxEi}  ,  {0,  Mauc[Sp]}},  Plot  Joined ->  True]  ; 


Plus@@Sp 

144.524 

■  Let's  fit  the  secondary  emission  data 

Lsd  =  ListPlot [Transpose [ {ExTest,  YxEMuData} ] , 

PlotRange ->  {{2,  6},  {1,7*10“*,  2,2  *10"®}},  DisplayFunction ->  Identity]  ; 
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YsMu  [x_]  =  Fit  [Tremsposo  [  {ExTest ,  YxEMuData}  ]  ,  }  /  x] 

2.65899  x  10"^ 

xi/4 

LsFit  =  ListPlot [Transpose [ {ExTest/  YsMu[ExTest] }] ,  Plot Joined  ->  True, 

PlotRzoige ->  {{2,  6},  {1.7*10*®,  2.2*10*®}},  Di splayFunct ion ->  Identity]  ; 
Show[{Lsd,  LsFit},  DisplayFunction ->  $DisplayF'unction]  ; 


Ys[Ex_]  :=  YsMu[Ex]  *Ex*  muAbs[Si,  Ex]  ; 

(Ys[ExTest]  +Yp[ExTest])  *  10^ 

(*  visually  check  the  fit  to  the  original  graphed  data,  ok  *) 

{13.5946,  6.72798,  5.79254,  5.10083,  4.48134} 

Ss  =  Spec  In  Ys  [Eout]  ; 

Plus@@Ss 

13.8085 


■  Now  smear  the  photoemission  in  energy 

■  First  the  primary  photoemission 

SpecPrim  =  TeU^le [0,  {nPts}]  ; 

For[i  =  1,  i  ^  nPts,  i++, 
newPE  =  Sp|[i]]  /  i  ; 

For[j  =  1,  j  ^  i,  j++. 

Spec  Prim  I  jj  +=  newPE 

]] 
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ListPlot [Transpose [ {Eout,  SpecPrim} ] /  Plot Range ->  {{0,  MaxEi},  {0,  Max [SpecPrim] } } ]  ; 


Plus  @@  SpecPrim 
144.524 

■  Add  the  secondary  photoemission  (all  in  the  bottom  bin) 

SpecTot  =  SpecPrim  ; 

SpecTotll]]  +=  Plus©@Ss  ; 

ListPlot  [Transpose  [{Eout#  SpecTot}]/  P  lot  Range ->  {{0/  MeixEi),  {0,  Mauc  [SpecTot]  }}  / 
Plot  Joined  ->  True]  ; 


